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Abstract 

We present and analyze the performance of constellation precoded beamforming. This multi-input multi- 
output transmission technique is based on the singular value decomposition of a channel matrix. In this work, 
the beamformer is precoded to improve its diversity performance. It was shown previously that while single 
beamforming achieves full diversity without channel coding, multiple beamforming results in diversity loss. In 



^ ' this paper, we show that a properly designed constellation precoder makes uncoded multiple beamforming achieve 

00 . 

ff^ , full diversity order We also show that partially precoded multiple beamforming gets better diversity order than 

■<!::j- ' multiple beamforming without constellation precoder if the subchannels to be precoded are properly chosen. We 



propose several criteria to design the constellation precoder. Simulation results match the analysis, and show that 
precoded multiple beamforming actually outperforms single beamforming without precoding at the same system 
data rate while achieving full diversity order. 
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<^ . I. Introduction 

When the perfect channel state infomiation is available at the transmitter, beamforming is employed to 
achieve spatial multiplexing and thereby increase the data rate, or to enhance the performance of a multi- 
input multi-output (MIMO) system HI. The beamforming vectors are designed in |[2l, [l3]| for various design 
criteria, and can be obtained by singular value decomposition (SVD), leading to a channel-diagonalizing 
structure optimum in minimizing the average bit error rate (BER) [3]. Uncoded single beamforming, 
which carries only one symbol at a time, was shown to achieve full diversity order of MN where M 
is the number of receive antennas and N is the number of transmit antennas JU, |[5l. However, uncoded 
multiple beamforming, which increases the throughput by sending multiple symbols at a time, loses full 
diversity order over flat fading channels HI, [HI. 
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It is known that an SVD subchannel with larger singular value provides larger diversity gain. During 
the simultaneous parallel transmission of the symbols on the diagonalized subchannels, the performance is 
dominated by the subchannel with the smallest singular value. To overcome the degradation of the diversity 
order of multiple beamforming, bit-interleaved coded multiple beamforming (BICMB) was proposed [[6]|, 
II3. This scheme interleaves the codewords through the multiple subchannels with different diversity order, 
resulting in better diversity order. BICMB can achieve the full diversity order offered by the channel as 
long as the code rate Rc and the number of subchannels used S satisfy the condition RcS < 1 [IHl, [|9l. In 
this paper, we present an uncoded single and multiple beamforming technique that achieves full diversity 
order. This technique employs the constellation precoding scheme [ fTOl , [[TT|. [[T2ll . |fT3l , O, which is 
used for space-time or space-frequency block codes to increase the system data rate without losing full 
diversity order. We show via analysis and simulations that fully precoded multiple beamforming achieves 
full diversity order even in the absence of any channel coding. For this purpose, we derive an upper 
bound for the pairwise error probability of the precoded beamforming system. Several criteria to design 
the precoding matrix are proposed. Simulation results show that fully precoded multiple beamforming 
with a properly designed precoding matrix outperforms single beamforming at the same system data rate. 

The rest of this paper is organized as follows. The description of precoded beamforming is given in 
Section [III Section UlI] calculates the upper bound of pairwise error probability for different schemes 
of precoded beamforming. In Section |IVl several criteria to design the precoding matrix are proposed. 
Simulation results supporting the analysis are shown in Section |Vl Finally, we end the paper with a 
conclusion in Section |VIl 

II. System Model 

The MIMO channel H G C*^^^ is assumed to be quasi-static, Rayleigh, and flat fading, and perfectly 
known to both the transmitter and the receiver. The beamforming vectors are determined by the SVD 
of the MIMO channel, i.e., H = UAV^ where U and V are unitary matrices, and A is a diagonal 
matrix whose s*'' diagonal element. As G M, is a singular value of H in decreasing order. When S 
symbols are transmitted at the same time, then the first S vectors of U and V are chosen to be used 
as beamforming matrices at the receiver and the transmitter, respectively. In Fig. \T\ which displays the 
structure of constellation precoded beamforming, U and V denote the beamforming matrices picked from 
U and V. Depending on S and the number of symbols precoded R, constellation precoded beamforming 
can be classified into three types as are described below. 



A. Precoded Single Beamforming 

In the Precoded Single Beamforming (PSB) scheme, a precoder 6^ generates a symbol y from an 
it! X 1 modulated symbol vector x, y — 0^x. We assume that each of the R members of x belongs to a 
normalized signal set x C C of size \x\ — 2"*, such as 2"* -QAM, where m is the number of input bits to 
the Gray encoder. Due to the beamforming vectors employed at both of the ends, the precoded symbol is 
transmitted on the subchannel with the largest singular value. Hence, the detected symbol r is written as 

r = Ai6>^x + n (1) 

where Ai is the largest singular value of the channel matrix H, and n is an additive white Gaussian noise 
with zero mean and variance A^o — N/SNR. H is complex Gaussian with zero mean and unit variance, 
and to make the received signal-to-noise ratio SNR, the total transmitted power is scaled as N. The 
maximum likelihood (ML) decoding of the detected symbol is given by 

X = min |r — Ai0^x|^. (2) 

The system data rate for PSB is r] = m - R bits/channel use. 

B. Precoded Multiple Beamforming 

In this scheme, > 1 modulated symbols are simultaneously transmitted on the subchannels with the 
largest S singular values. The S x 1 symbol vector x whose elements belong to x precoded by a 
square precoding matrix 0. The constellation precoder can be expressed as 



= P X 





Is-R 



(3) 



where @ is R x R constellation precoding matrix that precodes the first R modulated symbols of the 
vector X, and P is an 5" x 5" permutation matrix to define the mapping of the precoded and non-precoded 
symbols onto the predefined subchannels. When all of the S modulated symbols are precoded (R = S), we 
call the resulting system Fully Precoded Multiple Beamforming (FPMB), otherwise, we call it Partially 
Precoded Multiple Beamforming (PPMB). For FPMB, P can be replaced by the identity matrix, and 



= 0. The S X 1 detected symbol vector r at the receiver is written as 
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r = A50X + n (4) 

where is a diagonal matrix whose elements are the first S singular values of A, and n is an additive 
white Gaussian noise vector. The ML decoding of the detected symbol is given by 

1 1 2 

X = min ||r — As0x|| . (5) 
The system data rate for precoded multiple beamforming is r] = m ■ S bits/channel use. 

in. Diversity Analysis 

In this section, we will calculate the diversity order by analyzing the pairwise error probability (PEP) 
between the transmitted symbol x = [xi ... xs]^ and the detected symbol x = [£i ... xs]^- 

A. Precoded Single Beamforming 

For the ML decoding criterion of Q, the instantaneous PEP can be expressed as 

Pr (x ^ X I H) = Pr (|r - Ai6>^x|2 > |r - Ai^^xl^ | H) 

= Pr (/5 > |Ai6»^(x-x)n H) (6) 

where /3 = — Ai[0^(x — x)]*n — \in*9^{x. — x). Since /9 is a zero mean Gaussian random variable with 
variance 2A'oAf |0^(x — x)p, @ is rewritten as 



where Q{-) is the well-known Q function. By using the upper bound on the Q function Q{x) < ^e^^'^/^, 
the average PEP can be expressed as 



Pr (x ^ x) = ^ [Pr (x ^ X I H)] 



< E 



exp 



1 / A2|6l^(x-x^i2xn 



.2 V 4A/'o 



(8) 



Assume that \0 (x — x)p > for a distinct pair of x and x. Previously, in dH, [[8l, and BUl, we showed 
the closed form expression of ([8]). We provide a formal description of the result from [[8l, flU below. 



Theorem 1: Consider the largest S < min(M, A^) eigenvalues /i^ of the uncorrected central M x N 
Wishart matrix that are sorted in decreasing order, and a weight vector a = [ai ■ ■ ■ as]'^ with non- 
negative real elements. In the high signal-to-noise ratio regime, an upper bound for the expression 
£^[exp(— 7 ^^^^ agfis)] which is used in the diversity analysis of a number of MIMO systems is 

s 



E 



exp I -7 ^ asfi 



''S 



<C(«™„7)-^^-^^^^^^-'^^^ (9) 



where 7 is signal-to-noise ratio, C is a constant, amin = min{ai, ■ ■ ■ , as}, and 6 is the index to the first 
non-zero element in the weight vector. 

Proof: See m, [HI. □ 

Applying Theorem [T] to ([8]) where S = 1, 5 = 1, and amin = — x)P, we get the upper bound 
for the PEP as 

Pr (X - X) < i (i^^^^^^A^^) ■ (10) 

Therefore, it is easily found that PSB achieves full diversity order of MN once it satisfies the condition 

|6>^(x-x)|2 > (11) 

for any distinct pair of x and x. The method to design the precoding vector will be described in Section 
IIYI 

B. Fully Precoded Multiple Beamforming 

By using the same approach in PSB, we get the upper bound to the instantaneous PEP for precoded 
multiple beamforming as 

Pr (x ^ X I H) = Pr (||r - AsGxf > ||r - AcjGxf | H) 
1 / ||As0(x-x)f\ 

Let's define d = [di ■ ■ ■ dsY as a Euclidean vector that results from the precoded symbols of a distinct 
pair X and x. Then, d = 0(x — x), and the absolute value \di\ can be interpreted as a Euclidean distance 
between the symbols belonging to a new constellation transformed by the i*^ row vector of from the 



original constellation. For FPMB, the average PEP is expressed as 



Pr fx ^ x) < E 



2exp 



4iVo 



v 



(13) 



Applying Theorem [T] to (fT3l) . we get the upper bound to PEP as 

^{M-S+l){N-5+l) 



Pr (x ^ x) < C 



dr, 



AN 



-SNR 



(14) 



where C is a constant, dr, 



mm{\d 



2 

1| 5 



of the (squared) Euclidean distance vector [\di 



1^5^}, and S is an index to the first non-zero element 
■ ■ ■ Idsl"^]. Therefore, FPMB also achieves full diversity 

l^f (x — x)p > for any distinct 



order if 5 from any distinct pair is equal to 1, which implies that 
pair, where 6^ is the first row vector of 0. The way to build the precoding matrix will be described in 
Section HVl 



C. Partially Preceded Multiple Beamforming 

The partial precoding scheme divides the modulated symbols into two groups of symbols, i.e., precoded 
and non-precoded symbols. Through the permutation and the grouping, the numerator of the exponent 
term in (fT2l) can be represented as described below. For this purpose, let's define bp = [bp{l) ■ ■ ■ bp{R)] 
as a vector whose element bp{k) is the subchannel on which the precoded symbols are transmitted, and 
bp{k) < bp{l) for k < l.ln the same way, b„ = [6„(1) ■ • • bn{S — R)] is defined as an increasingly ordered 
vector whose element 6„(/c) is the subchannel which carries the non-precoded symbols. By reordering the 
resulting vector A50(x — x) for a simpler representation of A^P, the numerator of the exponent term 
in (fT2l) can be expressed as 















A5e(x-x)|| = 


A^P 








(x-x) 











-R 








' \ 






















A„ 







'^S-R 





Ap0(xp - Xp) 



(15) 
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where Ap and are the Rx R and [S — R) x [S — R) diagonal matrices whose elements are the ordered 
singular values corresponding to the subchannels of the vectors bp and b„, respectively, and similarly 
Xp = [xi ■ ■ ■ xr], x„ = [xn+i ■ ■ ■ xs], Xp = [xi ■ ■ ■ xr], x„ = [xr+i ■■■ xs]. By plugging ([B]) in 
we get an upper bound to PEP for PPMB as 



FtU^x) < E 



K 



^0 



(16) 



where 

R S-R 

= ^ '^bp{i)\di\'^ + ^ '^b„{i)\^R+i ~ ^fl+iP (17) 
i=l 1=1 

and (ij is the i^^ element of a Euclidean vector d = 0(xp — Xp). Let's assume that the constellation 
precoding matrix meets the condition of FPMB to achieve full diversity order. Since (fTTl ) has the 
closed form expression similar to (fT4l) as described in FPMB, 6 value needs to be obtained from a 
composite vector with such kind of elements as |c?jp and \xR+i — x/j+jp, to observe the diversity behavior 
of a given pairwise error. In addition, a different pair can lead to different diversity behavior. Therefore, 
we need to get the maximum 6 out of all the possible pairwise errors to decide the diversity order of a 
given PPMB system. 

All of the distinct pairs of x and x can be divided into three groups in terms of Xp, Xp, x„, and x„. 
The first group includes the pairs that have Xp = Xp, and the second group comprises the pairs satisfying 
Xp 7^ Xp, but x„ = x„. Finally, the last group is consisted of the pairs that Xp ^ Xp, and x„ 7^ x„. We will 
present the method to calculate 6 for a pair of each group, and to find 6max for each group. 

Since the vector d is a zero vector for the first group, the first summation of k in (fTTl ) is zero, resulting 
in 5 being equal to the minimum of b„. By considering all of the possible pairs, we can easily see that 
&n(l) < 5 < bniS — R). Therefore, the maximum value is 61 = bn{S — R) which corresponds to the pair 
satisfying Xj = Xj for all i except i = bn{S — R). For any pair in the second group, the term with the first 
singular value survives in k, according to the inherited property of the constellation precoding matrix, 
i.e., > 0. However, the second summation in k disappears since x„ = x„. Therefore, the maximum 
value of this group 82 = bp{l). Now, for the third group, both summations in k exist. Then, S is chosen 
to be the smaller value between the minimum of b„ and &p(l). In the same manner as was already given 
in the analysis of the first group, the maximum of the minimum of b„ can be found to be bn{S — R). 
Therefore, the maximum 5 for this group is ^3 = max{6p(l), bn{S — R)}. Finally, Smax can be decided 
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as 



5, 



max 



max{5i, 62, S3} 



max{6„(S' - R), bp{l), max{6p(l), 6„(5' - R)}} 



max{6p(l), bn{S - R)}. 



(18) 



Example: We provide the diversity analysis of the 4x4 partially precoded multiple beamforming system 
with 5 = 4 and R = 2. In this example, we assume that the precoded symbols are transmitted on the 
subchannel 1 and 3 while the non-precoded symbols are transmitted on the subchannel 2 and 4. Then, 
this configuration gives bp = [13], and b„ = [2 4]. By following the result in (fTSi) . 5max is equal to 
max{l, 4} = 4, leading to the diversity order of 1. The pairwise errors, satisfying xi = Xi, X2 = £2, x^ = 
X3, but X4 7^ £4, inflict loss on the diversity order of this system. Table IJ summarizes the diversity order 
analysis for all of the possible combinations of the 4x4 partially precoded multiple beamforming system. 



A. Precoded Single Beamforming 

The optimum precoding vector should maximize the array gain as well as meet the condition for 
achieving full diversity order. Based on (flOl) . the optimum precoding vector d^pt can be determined by 
solving the following equation: 



subject to the power constraint E\6^yi\'^ = E [|x||^ = R. The problem in (fT9l) coincides with that of IfTSl 
which addressed multi-user coding. In [[TSl . a number of users are assumed to transmit simultaneously 
at the same power, each of whom uses a rotated version of QAM symbols. Thus, the noiseless received 
symbol belongs to a new constellation containing the sum of each rotated QAM symbols. For a simple 
construction, the rotation vector is defined as 



IV. Precoder Design 



Sopt = arg max min 1 0^ (x — x) | 



(19) 



d 



1 e^<^ e^'^"^ e^'^"^ 




(20) 



where j = ^/^. The optimum rotation angles for R < 7 are found by searching that maximizes the 
minimum squared Euclidean distance of the new constellation. Hence, we can apply the result of [fTSll to 
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preceded single beamforming. 

B. Preceded Multiple Beamforming 

To design the preceding matrix, we establish various design criteria. Since we focus on even power 
distribution on each subchannel, the precoding matrix is restricted to be a unitary matrix which preserves 
the power. 

1) Maximization of the Minimum Euclidean Distance Among the S Coordinates, This criterion 
minimizes the upper bound to PEP of (fT4l) by maximizing dmin as 

@opt = arg max min dmin (2 1 ) 

subject to the power constraint i?||0x|p = ||x||^ = R. Since the analytical solution of this problem is 
unavailable, computer search can be used for small S and small constellation sizes. For this purpose, 
we employ the parameterization method of complex unitary matrices in tfTTTl . Li this method, any S x S 
unitary matrix S can be written as 

S = D H Gkii^pki^Pki) (22) 

l<k<S~l,k+l<l<S 

where D is an S* x S* diagonal unitary matrix, t/j^i E [—n,n], p^i E [—n/2,n/2], and Gki {fpkh pki) is 
a complex Givens matrix, which is the S x S identity matrix with the {k, ky^, (/, /)*^, {k, Z)*^, (I, kY^ 
elements substituted by costpki, costpki, e~-"''=' sIht/'h, and — e-'''''' sin^fc^, respectively. The optimization 
of the diagonal unitary matrix D is not necessary since dmin is the squared absolute value of an element 
which includes the diagonal entry of D with the magnitude equal to one. Therefore, we need to optimize 
only S{S — 1) parameters of {tpkh <Pki}- The optimum precoding matrix found by this method automatically 
satisfies the full diversity order condition since \di\^ > dmin > 0. 

2) Maximization of the Minimum Euclidean Distance of the First Coordinate, $2-' For large S, the 
squared Euclidean distance of i*^ coordinate is smaller than A^. Thus, the term A^|(iip takes the 
biggest portion of the summation in (fT3]) . This fact leads to an idea that maximization of the minimum 
Euclidean distance of the first coordinate |c/ip can lower the PEP. In this criterion, we solve the following 
equation as 

®opt = arg max min |^^(x — x)|^ (23) 

Vx^x 
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subject to the power constraint .EHQxlp = E = R. Since it is difficult to solve (|23l) in a tractable 
way, and the optimization equation is the same as (fT9l) . we propose a method that adopts the result of 
(fT9l ). In this method, we use the optimum precoding vector (fT9l) in &opt as the first row vector. To build 
a unitary matrix, we utilize the S'-point inverse fast Fourier transform (IFFT) matrix F5, whose {l,m) 
element is given by v^exp(j27r(/ — l)(m — The IFFT matrix provides two properties we can 

use; the elements of the first column vector are all ones, and the IFFT matrix is a unitary matrix. By 
constructing a precoding matrix as 

eopt = F|diag(^opi) (24) 

where Oopt is the precoding vector obtained from (fT9l ). we can see that the unitary matrix Gopt contains 
Oopt as the first row vector. This method guarantees full diversity order since > 0. 

3) Maximization of the Geometric Mean, $3.- Since the summation in (fT3l) consists of many terms, the 
previous optimizations which optimize only one term may not necessarily be the best criterion. Maximizing 
the arithmetic mean or the geometric mean of |c/jp values are attractive candidates which consider all of 
\di\^ values simultaneously. Between the two, we choose the geometric mean since the arithmetic mean 
does not necessarily guarantee the full diversity order condition. However, the geometric mean meets the 
condition on \di\^ > since the geometric mean of the optimum precoding matrix in this sense will be 
larger than zero. Therefore, the optimization is based on 

s 

©opt = arg max min TT 1 0j'(x — x)p/'^ (25) 

i=l 

subject to the same power constraint as the previous criteria, and Oj is the i*^ row vector of the precoding 
matrix 0. Authors in IfTTTl introduced an algebraic method to construct the precoding matrix in the space- 
time diversity system. According to [11], the unitary precoding matrix can be written as 

©opt = F|diag(l,a,--- ,cr^-i) (26) 

where a = e^'^'^^^, and the method to determine P is available in IfTTTl . 

V. Simulation Results 

To illustrate the analysis of the diversity order in Section [nil we now present simulation results over 
various channel dimensions. Fig. |2l [3l and |4] show bit error rate (BER) performance of SB, PSB, and 
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FPMB. The curves with legend FPMB $i, ^3 are generated by the preceding matrices based on 
each criterion in Section |Wl For a fair comparison between the different schemes, the system data rate 
1] is set to 4, 6, 8 bits/channel use. Throughout the figures, PSB and FPMB are shown to achieve full 
diversity order since the slopes are parallel to that of SB which is known to achieve full diversity order. 
A comparison between SB and PSB reveals that SB outperforms PSB for any channel dimension. The 
reason is that the array gain which can be observed by PEP is related to the minimum squared Euclidean 
distance, and the minimum squared Euclidean distance of the new constellation generated by the preceding 
vector is smaller than that of the SB constellation. For example, the minimum squared Euclidean distance 
of ?7 = 4 PSB normalized constellation is 0.27, while that of normalized 16-QAM SB is 0.4. 

Contrary to the case of 2 x 2, FPMB outperforms SB for larger channel dimension. In the case of 3 x 3, 
FPMB $3 gives 2 dB gain over SB. A bigger gain of 6 dB is observed by FPMB $3 for the case of 
4x4. Bigger gain for larger dimension can be explained by a comparison of the instantaneous PEPs of 
SB in BU and FPMB in (fT3l) . It can be stated that a larger number of singular values lead to a bigger 
array gain. We also find that the maximization of the geometric mean results in better performance than 
the others for larger channel dimension. 

Simulation results to support the diversity analysis of 4 x 4 5* = 4 PPMB in Table |I] are provided in 
Fig. [51 We find that the simulation results follow the diversity orders in Table HI The curves with the 
same diversity order give different array gain depending on the subchannel selection of preceded symbol 
transmission. BER at high SNR for bp = [1 4], [2 4], [1 2 4] are the same, leading to around 10 dB larger 
gain than that of bp = [3 4]. This can be explained by the fact that PPMB with the smaller number of 
pairs (causing the worst diversity order) provide larger array gain. Since the subchannel 3 transmitting the 
non-precoded symbol dominates the performance loss for bp = [14], [2 4], [12 4], only one pair, satisfying 
X3 = X3 and Xi ^ Xi for i = 1, 2, 4, is related to the worst diversity order. On the other hand, for bp = [3 4] 
where the subchannel 3 also dominates the performance loss, half of the total possible pairs, such that 
xi = xi, and X2 = £2, and X3 7^ £3, or X4 7^ x^, make the worst diversity order. This fact also applies to 
the case of diversity order of 9. 

VI. Conclusion 

In this paper, we introduced a preceded beamforming system that achieves full diversity order. For 
the analysis, we calculated an upper bound to the pairwise error probability, assuming that the receiver 
decodes the transmitted symbols based on maximum likelihood decoding. We established several criteria 
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to design the precoder. We also provided simulation results that support the diversity analysis of precoded 
beamforming. We showed in particular fully precoded multiple beamforming with a proper precoding 
matrix outperforms single beamforming at the same system data rate while achieving full diversity order. 
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TABLE I 

Diversity order of 4 x 4, S = 4 partially precoded multiple beamforming system 
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Fig. 2. BER vs. SNR comparison for 2 x 2 16-QAM SB, 4-QAM R = 2 PSB, and 4-QAM S 



= 2 FPMB. 
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Fig. 3. BER vs. SNR comparison for 3 x 3 64-QAM SB, 4-QAM R = 3 PSB, and 4-QAM 5 = 3 FPMB. 
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Fig. 4. BER vs. SNR comparison for 4 x 4 256-QAM SB, 4-QAM R = A PSB, and 4-QAM 5 = 4 FPMB. 
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